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Abstract 

Perturbation of FRW spacetime is carried out in NP formalism. 
The equation governing the scalar, vector and tensor modes take on a 
very simple and transparent form. All of them can be combined in one 
master equation for all helicities. The solutions for the closed, flat and 
open FRW are analytic continuation of the same functions, so only 
the solutions in the closed model are described. The scalar equation 
is the same as that of the conformally coupled massless Klein-Gordon 
field, the vectorial ones are the same as Maxwell equations, and the 
tensorial ones are for spin-2 fields. The corresponding eigen-functions 
are all determined, and in particular, the Green's function and the 
Lienard-Wiechert type potential also solved. These solutions reduce 
to the familiar form in fiat space. 

1 Introduction 

Motivated by the success of the Newman-Penrose(NP) method in investi- 
gating the perturbation of various black hole geometries, we thought that 
the application of this method to study the perturbation of the Friedmann- 
Robertson- Walker spacetime may give some further insight into structure 
formation. One of the authors [H [2] had used this formalism to study free 
Maxwell and Dirac fields in FRW space. After setting up the notations be- 
low, we write down the perturbation equations in the next section and solve 
them. Section 3 discusses the scalar perturbation, and the Green's function 
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as well as the Lienard-Wiechert type potential in closed space are found In 
the final section we discuss some consequences. 

Writing the line element as ds'^ = \dif — dr'^ — 5^ [dO'^ + sin^ 6d(f)'^)] 

with S = — ^ =, we may devise the null tetrad to be = [1, — 1, 0, 0], = 
^[1, 1, 0, 0], = —^[0, 0, 1, i sin6']andthecomplexconjugate m^, to express 
the directional derivatives as D = l^d^ = A = n^d^ = — S = 

= 7^^"'^* = = vts^"- ^ith = J, T 1^ and £± = 

W "F iil^e^- Also, the non- vanishing spin coefficients[3] are given by /3 = 
-a = 7 = 2n = ^ - ^, -a^p = ^ + where the overdot and 

prime denote derivatives with the conformal time r] and r respectively. We 
will also have to use the operators T>^ = T>^ + and = + ncotO. 

With the change of variables ^/Kt] = x, and sin ^V^'"^ = cosec^, we get 

sin (^y/Kr^ T>^ — )■ iA^„ = i ± ^j^^ — ricot,^^, and A looks the same as 
C. 

Of the NP quantities representing the curvature of the spacetime, all 
the Weyl scalers (\E''s) vanish in the homogeneous and isotropic FRW back- 
ground, and the non-vanishing background Ricci scalers are given through 
R = 6 (^) , *oo = i^n = ^$22 = ^ (Ka^ + 2d^ - aa) = g(p + p). In 
the work below, we just consider the closed Universe with i^' = 1, as the 
other two, K = (0,-1) turn out to be analytic continuation of the same 
solution. 



2 Perturbation of the space-time 

Now we look at the perturbation equations. In combination with the \E''s 
that vanish in the background, we should take the unperturbed directional 
derivatives and spin coefficients for first order perturbation equations. Then 
we can write the first order perturbed Bianchi identities as 

tAf^f± C^_,^f_,= Qf, (1) 
tAts<^f_,T Ct<l>t= Uf (2) 

where either the upper or lower signs are taken once at a time, the Q' s and U' s 
are the source terms comprising the perturbed energy-momentum; the <l>'s 
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the gravitational or tensorial perturbations, $f 



are related to the Weyl scalars \1'' s, with $^ 




3, 



representing 



+ 



repre- 



senting the vectorial or vortical perturbations, and $o = the scalar 

perturbations. The equations for the vectorial perturbations are the same as 
the Maxwell equations discussed in Ref.[lJ. It is an interesting proposition 
that the vortical perturbations arise from the electrodynamic interaction be- 
tween the Maxwell and Dirac fields. Simlarly, the equation governing the 
scalar perturbation is the same as that of the conformally coupled massless 
Klein-Gordon field. There are 10 degrees of gauge freedom in the problem as 
we have formulated, four relating to the general co-variance of general rela- 
tivity and six corresponding to rotations of the tetrad frame as discussed in 
Ref. [3J. These conditions can best be used to simplify the source term, but 
we will not be concerned with the sources in this work, and will only discuss 
the eigen- modes of these perturbations. Inspecting Eqs. ([T]and|2]), we see 
that all the perturbations are related to each other, in that one contributes 
to the other; existance of any one single type of perturbation gives rise to all. 

Eqs, (P and |2]) can be decoupled by operating on the former with i A^_^ 
and the latter with ^Cf and adding the two. The result is 



Thus the variables are found to be separable. It is easy to establish the 
operational identity = C^^^CZs — ^^s, and a similar one for A. Letting 

p = ±s for the appropriate helicities, all the perturbation equations contained 
in Eq. ([3]) can be written in a master equation as as 



Changing the sign of p together with the signs on the operator gives the 
equation for fields with helicity reversed. The eigenf unctions are separable, 
and we may write$p =p Z XpY where pY are the angular part represented 
by the the spin-weighted spherical harmonics satisfying the 



ct^pCZppvr = -{i-p){i+p+i)pYr; 

the solution[l] is pY,"* {6, 0) = Are™<^(l - cos^)'^(l + cose)"^ P^^:'^''"^ , 



allowing us to identify the spin-weighted spherical harmonics as the spher- 
ical harmonics formed with the Jacobi polynomial Pn"'^\ It can easily 



(3) 



[-^t+p^-p + ^t+p'C-Zp] <^p = Jp. 



(4) 
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be established that pYj""* = (-l)"^+P_pY;-™, and oY{^ = Y^^ are just the 
usual spherical harmonics. These are normalized to give Yjp- *pl^™(idf2 = 

Sii'Smm', and they form a complete set in that Ylim p^™(^')*p^z™'(^)~^(^" 
Koornwinder[5j developed an addition formula for Jacobi polynomials, 
and there is a generalized addition theorem for pY/^ given in Ref.|6j, but 
we can devise a degenerate form for our purpose by noting that 5 (fi — Q') 
can depend only on the angle /3 between the directions Q and Q' where 
cos (3=cos e cos ^'+sin 6 sin 6' cos (0 - 0')- So, multiplying 5(fi-(^')=E« bi pY'"^ 0) 
on both sides by pYj7^ 0) and integrating over the solid angle to find 
bi =p Yi~^ (0, 0), we may write 

J2 m m = pY^* (o, o) pYr (/?, o) . (5) 

m 

For the eigen- functions of A we may substitute pZ'j! [rj, r) = e~^'^^ pR^ (r) 
to separate the temporal part and write the operator = ± + racot,^ 
giving the the radial eigen-equation 

At^pAZp pRt = {k-p){k + p+l) pRt (6) 

that determines the eigen-functions to be 

pRt = Ar-^-' (1 - cosO""^ (1 + cosO""^ PtV-r^"^ (cosO 

= Az"-'=-i+^' sin'^ r e-^P^P^XP'-"+^') {% cot r) ; (7) 

A is the normalization constant, and the parameters of the Jacobi polyno- 
mials are chosen to make the function regular at r = and tt. Ref. [7] has 
proved the non-Hermitian orthogonality of Jacobi polynomials with general 
parameters. In our case, we find that fiJ^°^ ip^k'p^k) (cos^) = 5kk' gives 



us A 



iw— 1 



(2fc+l)(a;-fc-l)!(A:+p)!(A:-p)! 
7r(tj+fc)! 



So pi?^* = Both the angular 



and radial functions satisfy the spin raising and lowering operations 



^±ppYr' — ^L±p(^p^i)Yl^ a,nd 

^±ppZk = ±-^±p(pTi)^fc) (8) 

where the eigenvalues are Lp = ^^/ {I — p) {I + p + 1) and Kp = ^y {k — p) {k + p + 1). 
In Eqs. ([8]), the upper signs lower the helicities while the lower sign raise 
it. Also the eigen-values satisfy L_p{oYK_p) = Lp^i (oiKp^i). Some of the 
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eigen-functions representing the various modes of perturbation are displayed 
in fthe accompanying figures. The scalar perturbations, ([T]), represent the 
density perturbations. These give rise to the structures. Fig. (|2]) show the 
vectorial perturbations. Thes are responsible for generating rotational or 
vortical effects on the perturbations. The next mode shown in Fig. ([3]) are 
the tensorial perturbation representing the gravitational radiation content. 
All these modes are inter-related; e.g., any scalar dennsity enhancement due 
to gravitational coUaps generates a rotational component and induces the 
release of binding energy through gravitational radiation. 



3 Scalar Perturbation 

Let us look at the scalar perturbation in some greater detail. To solve the 
inhomogeneous equation with source, we can work out the Green's function 
by writing Eq. (jl]) for p = with the delta function source: 

O O -I "I -| 

G = -^6 [7] - i) 5{r- r') 6 {9 - 9') 6{<p-<p') 
sin r 

(9) 

The temporal eigen-functions are e the angular are the spherical har- 



-|- .2 ''-1 ''-0 



Qj,2 Q^2 gjj^2 ^ 

he temporal eigen-1 
monies o^;™ = and the radial ones are the normalized and appropriately 



weighted Gegenbauer polynomial \fN sin'"^^ r C^^^. Now we can immediately 
write down the eigen-function expansion of the Green's function as[8] 

1 /"°° du] 



k,Lm ' 



(k + l + l) 

xN (sinr sinr') i+^Cl+^ (cosr) C^+i (cosr') (10) 

Firstly, the integral over u can be done by the method of residues to give 
IZo = [uj{k + l + l){T]- r]')] Next, the addition 

theorem for spherical harmonics can be used to write Y[^* {Q!) [Vt) = 
2|ti(^^(i/2) (^cos /3) where cos/3 is the angle between the directions and fi'. 
Then we are left with 

^ oo 
fc=0 

^ 2^'(2/ + l) v.' jk-iy. , . . ,^(m), , 

X (ic _\- ] 4- ^'\\ (smr smr j C^_; (cos r) G^_, (cos r 

1=0 -I- -I- 
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Here, we may use the addition theorem for Gegenbauer polynomials |¥] to 
COS r cos r' + sin r sin r' cos /5 . Thus, 



replace the second summation with Cn (cosp) = s"pp ? where cosp 



oo 

G = — {cos k [p - [t] - 7]')] - COS k[p-{ri- 7]')]} 

47r^ sm p ^ — ^ 
^ fc=i 

1 {5[p-(r^-V)]-5[p+(r/-r/')]}- (H) 



47r^ sin p 

The first term within the braces is the retarded, and the second is ad- 
vanced. Green's function. In the limit to the fiat case (r, r') — )■ 0, we find 

sinp — )■ p, and cosp — )■ 1 — ^ = ^f' cos P ^ p^ — )■ 
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r — r ; hence we recover the familiar Green's function from classical elec- 
trodynamics. With Eq. (ITT]) , we can consider the potential s [r]' , r' , 6' , (f)') = 
6[r' — ( [f]')] 6 [6' — [f]')] 5 [0' — 0c iv')] of ^ point scalar perturbation mov- 
ing along the trajectory {( {i]) ,6(;{ri) ,(j)t^{'r])) to find a retarded Lienard- 

Wiechert type potential. Thus, $o {p,v) = I Gs dr'd9'd(j)'dr]' = f H iP (V) - - ^')] 
The integral over T^'can now be easily done to give 



47r sin [p(r70] 1 + 



(12) 

v'=v-pW) 



where the variables are to be evaluated at the retarded time rj' = t] — 
p {rj')] here, cos [p {j]')] = cosr cos [r' {r]')] + sinr sin [r' (77')] cos [f3 (77')] and 
cos [(3 {rj')] = cos 9 cos [9' (r7')]+sin 9 sin [9' {r]')] cos [0 — (V)]) with the primed 
co-ordinates being the location of the source that produce the potential at the 
un-primed location of the observer; also the right hand side is to be evaluated 
at the retarded time. Without loss of generality, we may simplify the notation 
by locating the observer at the origin, in which case cos [p {1]')] = cos [r' [r]')]. 
Again it can be easily checked that Eq. (fT2|) reduces to the familiar form in 
the limit r — )■ small. 

4 Conclusions 

In this work we have written the equations of perturbation of the FRW space- 
time in Newman-Penrose formalism. We find that the tensorial, vectorial and 
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scalar perturbations decompose into forms that reveal the spin content spin 
content very transparently. The perturbation modes are also inter-related, in 
that the presence of any one mode gives rise to the rest. The eigenfunctions 
of all these modes of perturbation can be solved in terms of appropriate spin- 
weighted spherical harmonics and radial functions of similar spin content. 

The scalar mode gives rise to the density perturbations that are respon- 
sible for the large scale structures. We are able to solve for the Green's 
function of this mode, and the result reduces to the familiar form from clas- 
sical electrodynamics in the limit to fiat space. With the Green's function, 
we should be able to solve for the potential driving the scalar perturbations. 
Also, the scalar perturbation \E'o = ^3 ^"3 ^ may be interpreted to represent 
the Newtonian potential, as in the Schwarzscild spacetime "i/o = —^. 

The vectorial perturbation equations are the same as Maxwell equations. 
So it is worthwhile to investigate the possibility that the vortical modes and 
the rotational effects on the large, and even small, scale structures could have 
their origin in electrodynamic interactions in very early Universe. Although 
the free Dirac and Maxwell equations have been solved in another work, it 
is also necessary to to incorporate the quantum corrections into the model. 
All the known massive particles are fermionic, so study of the Dirac current 
and energy-momentum as well as their interaction with the Maxwell field will 
reveal many features of structures. These and similar problems will be taken 
up in future work. 

Regarding the tensorial perturbations, all that can be said at this mo- 
ment is that they must represent the gravitational radiation content of the 
Universe, and the gravitational waveforms radiated by the structures present 
in the Universe. 

The perturbation equations, Eq. (fT][2|) . show that source of any one mode 
of perturbation produces the other modes. So the whole picture will become 
consistent only in an integrated view. This work can be considered only as 
an initial phase. 
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Figure 1: The scalar radial eigen-f unctions. The dashed curves are the nor- 
malized functions and the solid curves are the squares of them. The upper 
two are n — u — k — 1 — 0, the lowest possible state on the with k — on 
the left, and k — 6 on the right. The lower ones are for n — 6, with k — 3 on 
the left and /c = 8 on the right. It is seen that the higer k states are confined 
towards the middle, nearer to r — n /2. 
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Figure 2: The vectorial radial eigenfunctions. The dashed curves are the 
real part, the dot-dashed are imaginary and the solid are the square of the 
modulus, n = for the upper two graphs, with the lowest state of /c = 1 on 
the left, and A; = 4 on the right. The lower graphs are for n — 3, with k — 1 
on the left and A; = 6 on the right. Again, the higher /c's are confined nearer 
to the middle. 
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Figure 3: Same as Fig. ([2]), with same values of n for the wave forms of 
the tensorial gravitational radiation. The upper left is the lowest state with 
k = 2, and the right is = 6. The lower are k = 4 and 10 respectively. 
Higher k values show similar tendewncise as the scalar and vectorial. 
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